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Abstract 



We propose a plasma model for spectral statistics displaying level repul- 
sion without long-range spectral rigidity, i.e. statistics intermediate between 
random matrix and Poisson statistics similar to the ones found numerically 
at the critical point of the Anderson metal-insulator transition in disordered 
systems and in certain dynamical systems. The model emerges from Dysons 
one-dimensional gas corresponding to the eigenvalue distribution of the clas- 
sical random matrix ensembles by restricting the logarithmic pair interac- 
tion to a finite number k of nearest neighbors. We calculate analytically 
the spacing distributions and the two- level statistics. In particular we show 
that the number variance has the asymptotic form T?{L) ~ xL for large L 
and the nearest-neighbor distribution decreases exponentially when s ^ oo, 
P{s) ~ exp(— As) with A = 1/x = A;/3-|- 1, where /3 is the inverse temperature 
of the gas (/? =1, 2 and 4 for the orthogonal, unitary and symplectic symme- 
try class respectively). In the simplest case of /c = /3 = 1, the model leads to 
the so-called Semi-Poisson statistics characterized by particular simple cor- 
relation functions e.g. P{s) = 4sexp(— 2s). Furthermore we investigate the 
spectral statistics of several pseudointegrable quantum billiards numerically 
and compare them to the Semi-Poisson statistics. 
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I. INTRODUCTION 



Since the advent of random matrix theory there has been considerable interest in the 
statistical analysis of spectra [0-3. Two diametrically opposed statistical distributions have 
been found to be of universal relevance: the Poisson distribution, i.e. completely uncorre- 
lated levels, and the Wigner-Dyson distributions of random matrix theory. The prominent 
features of these distributions are conveniently characterized with the help of spectral ob- 
servables such as the nearest-neighbor spacing distribution P{s) and the number variance 
S^(L) (the variance of the number of levels in an energy interval of length L in the unfolded 
spectrum). The former stresses the correlations on a short scale, while the latter measures 
the stiffness of the spectrum, i.e. long-range spectral correlations. For the standard random 
matrix ensembles (orthogonal, unitary and symplectic symmetry labeled by f3 =1, 2 and 4 
respectively) the spacing distribution is approximately given by the Wigner surmise and the 
number variance increases only logarithmically with large L, 

P{s) = aps^ exp(— C/js^) 
S2(L)~-^log(L) (L^oo) (1) 

{ap and Cjs are determined by normalization and by requiring the mean level spacing to be 
one). Thus the Wigner-Dyson distributions are marked by level repulsion and long-range 
spectral rigidity. In contrast, for the Poisson distribution one has neither level repulsion nor 
spectral rigidity, 

P(s) = exp(— s) 
T?{L) = L. (2) 

In the present article we devise a discrete set of statistical distributions with properties 
intermediate between those of the Poisson and Wigner-Dyson distributions, namely level 
repulsion 

P{s) r^s" {s-^ 0) (3) 

paired with an exponential decay of the nearest-neighbor spacing distribution, 

P{s) ~ exp(-As) (s -> cx)) , (4) 

and a linear asymptotic of the number variance, 

S2(L) r^xL (L^oc). (5) 

Our distributions arise as the Gibbs distribution of a one- dimensional gas of interacting 
classical particles. It was noted already by Dyson that the distribution of eigenvalues Ej of 
random matrices is equivalent to the distribution of particles in a one-dimensional gas with 
inverse temperature (3 and the repulsive two-body interaction potential V{x) = — log|a;|, 
since the joint probability distribution can be written in the form 
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where the normahzation constant Z^q is the canonical partition function of the gas and 



W{E^, ...,Em) = Y: U{E,) + Vm - E,\) (7) 

i i<j 

is the total potential energy. The one-body potential U serves to confine the Ej to a finite 
stretch of the ii^-axis, e.g. U{E) oc E'^ for the Gaussian ensembles. It determines the mean 
particle (level) density, but it is irrelevant for the statistical correlations on the scale of the 
mean level spacing (random matrix universality) as long as U is steep enough to actually 
confine the levels 

Our plasma model of intermediate spectral statistics is defined by the same eqs. @ and 
(0), except that we restrict the second sum on the right hand side of eq. (|^) to a finite 
number k of nearest neighbors. It is this screening of the Coulomb interaction between the 
levels, which is the essential ingredient leading to the properties eqs. (^-(^• 

Our first motivation stems from the physics of disordered electronic systems. A 3- 
dimensional disordered sample undergoes a phase transition between an insulating and a 
metallic phase as a function of the disorder strength (Anderson metal- insulator transition). 
In the insulating phase the electron eigenfunctions are localized and since non-overlapping 
eigenfunctions are uncorrelated, the eigenenergies are Poisson distributed. In contrast, in the 
metallic phase the eigenfunctions extend homogeneously over the whole sample and overlap 
strongly which leads to a Wigner-Dyson distribution of the energy levels. Exactly at the 
transition point the electron eigenfunctions are extended, but strongly inhomogeneous (mul- 
tifractal). This leads to intermediate spectral statistics, which are believed to be universal, 
i.e. independent of the microscopic details of the disordered system 0. 

Fyodorov and Mirlin 0] calculated the overlap of two critical eigenstates with slightly 
different energy and found that the overlap is one if the energy separation, s, between the 
eigenfunctions is of the order of mean level spacing, while for larger s it decays as a certain 
power of s. In contrast, for Wigner-Dyson statistics the overlap is one for all values of s. 
Hence one may conclude that at the critical point of the Anderson metal-insulator transition 
only eigenvalues which are separated by at most a few level spacings interact strongly. 

A second motivation which led us to consider the short-range plasma model is the follow- 
ing. In the context of quantum chaos the statistical analysis of spectra plays a major role. 
In the semiclassical limit the quantum energy spectra of systems with integrable classical 
dynamics generally display Poisson statistics |Q, while the Wigner-Dyson distribution is 
associated with fully chaotic classical dynamics |^ . But there are systems which are neither 
integrable nor chaotic. In fTH] we examined numerically a few such systems and found that 



their spectral statistics exhibit all properties (||) - (|5|) typical for intermediate statistics. In 
particular, the simple expression 

P{s)=4se-^' (8) 

is an excellent fit to the spacing distribution of a subclass of systems considered. We further 
noticed that the short-range plasma model with nearest-neighbor interaction {k = 1) leads 
precisely to this distribution (and we checked that other correlation functions are also well 
described by this model). 

Another interesting class of plasma models with screened Coulomb interaction is the 
Gaudin model defined by eqs. (^ and (0) with the two-body interaction potential 
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V{x) = -]-\og-^^. (9) 

When a — >■ oo it reduces to the random matrix models, while in the limit a — it yields 
Poisson statistics. For (3 = 2 the model is solvable and all correlation functions can be 
written in a closed form [11,12]. They obey all three conditions - (j^) characteristic for 
intermediate statistics with the following constants (for [3 = 2) 

X=-(l-e-'^), (10) 
a 

and 

A=ir-Mi±^)*. (11) 

Here a = 27iap and p is the mean density of the levels. 

In general, in any plasma model with short-range interaction the nearest-neighbor spacing 
distribution will decay exponentially at large distances. It is the otherwise unusual property 
of long-range interaction in the standard random matrix models (which manifests itself in 
the summation over all pairs of particles in eq. (|^)) that is responsible for the exp(— cs^) 
decay of P(s) at large s. 

The main drawback of the Gaudin model is that no analytical solution is known for 
(3^2. In contrast, our model can be solved for arbitrary (3 (and arbitrary potential). 

We stress that our one- dimensional gas model is meant as a toy model which deserves 
interest, because it provides a natural discrete interpolation between Poisson and Wigner- 
Dyson statistics, it is analytically solvable, and it is physically motivated. It does not pretend 
to furnish an exact description of the critical statistics of a physical model (e.g. the Anderson 
model in the MIT point). Despite significant theoretical progress (see 0, [l^, ||2^- |^ and 



references therein) , the precise form of the critical distribution is not yet known though there 
exist arguments ( |^) that the main difference between the standard random matrix 
ensembles and the intermediate statistics is in the form of the correlation kernel, K{x,y), 
which determines the correlation functions in the random matrix theory 0. In the standard 
ensembles 

, sin7r(x — w) , , 

K{x,y) = ^ f. 12 

n{x-y) 

For intermediate statistics it has been argue that 

asin7r(x — -u) 

K{x,y)^ \ ^, 13 

vr smh a{x — y) 

where a is a parameter. When a — the standard random ensembles are recoved. Non-zero 
value of a characterizes the intermediate statistics. The detailed discussion of this type of 
critical behavior is beyond the scope of this paper. 

The paper is organized as follows. In the next section we present the calculation of the 
spacing distributions and the number variance in the one-dimensional plasma, first for k = 1, 
then for k = 2, and finally for arbitrary k. We compare the resulting distributions with our 
numerical results for the pseudointegrable billiards in section |T|. Finally, in the discussion 



of section IV we compare the short-range plasma model to other existing interpolations 



between Poisson and Wigner-Dyson statistics. 
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II. THE MODEL 



Instead of a linear one-dimensional gas as in eqs. @ and (^), where the levels are confined 
by a one-body potential U{E), we choose a circular geometry. This is convenient, since on 
a circle the levels are automatically confined and it is unnecessary to introduce an external 
potential. The mean level density is then constant and unfolding is trivial. At the same 
time the correlations in the unfolded level sequence are the same as in the linear geometry 
(in the limit of a large number of levels), just as the Gaussian ensembles of random matrix 
theory are locally equivalent to Dyson's circular ensembles. The method of calculation that 
we apply in this section goes back to Giirsey ||13| for the case of nearest-neighbor interaction. 



and to van Hove for the general case (see also ref. |jT5[). 

We consider particles (representing energy levels) on a circle of circumference L. We 
denote the positive spacings between neighboring particles by with j = 1 . . . N (see Fig. 
Hence, 

N 

For convenience we use a periodic index, i.e. Cj+N '■= ^j- We introduce an interaction 
between the particles via a repulsive two-body potential V{^) (eventually we will choose 
V{^) = — log 1^1 as in random matrix theory), but we let each particle interact only with its 
k nearest neighbors to the left and to the right {k = 1, 2, 3, . . .). 
The canonical partition function of this one- dimensional gas is 

N 
i=l 



Z^{L, /5) = 1 d^ . . . 1 d^^ 5(l - ^ 

*=i 

N 

exp(-/55:iy(0,...,0+fe_i)) , (14) 



where 



w{^„ . . . , 0+fc-i) = ^(0) + ^(0 + O+i) + 

• + + + (15) 

W{^j, . . . ,C,jjf-k-i) includes the interaction energy of the particle which is located at the 
left of with its k nearest right neighbors only to avoid double counting. Since for the 
application to spectral statistics j3 is not a free variable, but takes only the values 1, 2, and 
4, we suppress the /3-dependence of the partition function in the following. For brevity of 
notation we define 

/(0:=exp(-/3F(0) , (16) 
which for the case of ^(^) = In |^| becomes 

fio = (17) 

The joint probability distribution of n consecutive spacings then takes the form 
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P{^1, • • • , ^n) = -^^^ J ^^n+l ■■■J d^N S(^L 





N 

X n /(^.)/(^. + M ■ • ■ /(e. + • • • + w-i) • (18) 

In the following the variable s denotes distances measured in units of the mean spacing 
A = L/N. Our aim is to calculate correlation functions such as the n-th nearest-neighbor 
spacing distribution P{n, s) (the distribution of the distance s between two particles that 
have n particles between them) and the two-point correlation function R2{s) (the probability 
of finding any two particles at distance s). The latter is related to the former by summation 
over n, 

oo 

R2{s) = Y.P{n,s). (19) 

n=0 

Using eq. (|T8|) the n-th nearest-neighbor spacing distribution can be expressed as 

oo oo „ I -I 

P(n, s)=Jd^,...J dUi - A E X 



p{^i,...,Ui) ■ (20) 

The number variance can be obtained from the two-point correlation function with help of 
the relation 

L 

E'^{L) =L-2j ds{L-s){l-R2{s)) . (21) 



A convenient way to calculate the asymptotic behavior of the number variance for large L is 
to consider the asymptotic expansion of the Laplace transform of the two-point correlation 
function, 

oo 

g^it) =Jds R2{s) e-'' , (22) 



for small t. If 

g2{t) = j + ao + ait + 0{t^) (t ^ 0) (23) 



then we get from eq. (^TD 

j:'iL) = xL + C + OiL-') (L->oo), (24) 

where x = 1 + 2q;o and C = 2ai (for the determination of the constant term we have assumed 
that R2{s) falls off faster than which is true for the cases considered in this article). 
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A. Nearest-neighbor interaction (k = 1) 



We begin with the simplest case, where the interaction is restricted to nearest neighbors, 
so that the expression for the partition function eq. (|1^ simphfies to 



N N 



Zn{l) = I d^i . . . I de^v - e n fii 



(25) 



i=l j=l 

By Laplace transformation with respect to L we obtain 



dx /(x) e' 



-tx 



N 



git) 



N 



(26) 



The large N limit of the partition function can now be calculated by performing the Laplace 
inversion in saddle point approximation. We have 



1 

27Ti 

1 
27r^ 



C+JOO 



dt gN{t) e 



Lt 



C— lOO 

c+ioo 



gAr(tA+logg(t)) 



C— lOO 



Lc 



(27) 



where A = L/N is the mean level spacing and c is determined from the saddle point equation 

A + ^ = . (28) 

9[c) 

The expression for the joint probability distribution of n consecutive spacings eq. ([18[) 
reduces in the present case to 

. . . , e.) = ^^-^^;F^^^^^ n m) ■ 



Zn{L) 

Using eq. (pTl) and assuming n ^ N we obtain 



M6,...,en) = n -7;Te-^^^/(e,). 

j=i y\^) 



(29) 



Note that the factor g[c) ^ assures the proper normalization, 

oo 

j d^n • • • , -^n) = P(6, • • • , in-l) ■ 
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With help of eqs. ( pOD and ( ]29|) we find for the Laplace transformation of the n-th nearest- 
neighbor spacing distribution 



g(n,t) := J dse~^'P{n,s) 



g{c + t/A) 



-\ n+l 



9{c) 



(30) 



Using the relation (|T^ we calculate the Laplace transform of the two-point correlation 
function, eq. (p2|), 



The small-t asymptotic behavior is 



gjc + t/A) 

g{c)-g{c + t/A) ' 



(31) 



t 2A g'[c) 

, ?>g"{cf-2g\c)g^^\c) t 



12 g>{cy 



A2 



(32) 



where we have employed the saddle point equation (^). 

At this point we specialize to the random matrix interaction potential, i.e. we substitute 
eq. (p!?!). Its Laplace transform reads 



9{t) = 

and the saddle point equation (|^) yields 



r(/3 + i) 



[3 + 1 



(33) 



(34) 



By Laplace inversion of eq. (^) we then obtain the n-th nearest-neighbor spacing distribu- 
tion 



Pin^s) 



+ ^> «/3+n(/3+l) ^-{P+l)s 



r((/? + l)(n + l)) 

In particular, the nearest-neighbor spacing distribution for /? = 1,2,4 is 

P{s) = Ase-^' (/5 = 1), 

P(.) = |.^e-^^ (/3 = 2), 



(35) 



Pis) 



3125 
~2A' 



(/9 = 4) . 



The Laplace transform of the two-point correlation function eq. (0) becomes 

1 



(36) 



(37) 



8 



from which R2{s) can be computed by Laplace inversion. 

i?2(s) = 1 - exp(-4s) {(3 = 1) 

i?2(s) = 1 - (cos(3y3s/2) + 



V3 sin(3y3s/2))e-9^/2 ^^^2) 

4 



fc=0 



With help of eqs. ( P5D and and the small t expansion of eq. (^) we find the large L 



behavior of the number variance ~ xL + C with 

1 ^ /^(/^ + 2) 

^ + 6(l + /3)2- 

The spacing distribution and the two-point correlation function of the simplest model with 
k = j3 = 1 are very close to the corresponding distributions of a number of different dy- 
namical systems, whose spectral statistics can be calculated only numerically (see sects.pl 
and |rv| ). In lack of a deeper understanding of the spectral statistics of these systems, the 
plasma model is valuable in that it provides simple fitting distributions to these numerical 
results. We propose to call the statistics derived from the plasma model with k = (3 = 1 the 
Semi-Poisson statistics. This name was originally motivated from the fact that if one takes 
an ordered Poisson distributed sequence the nearest-neighbor spacing distribution of 

the new sequence {?/„} with ?/„ = (x„ + x„+i)/2 coincides with eq. (|) (of course the other 
correlation functions of the sequence {Vn} are different from those of the plasma model). 

It is interesting to note that again starting from the sequence and dropping every 
second level one obtains a new sequence with precisely the same statistical distribution as 
the plasma model with k = (3 = 1 More generally, retaining only every {p + 1)— th level 



of the sequence leads to the same statistical distribution as the plasma model with 

A; = 1 and/3 =p [0. 

B. Interaction between nearest and next-to- nearest neighbors (k = 2) 

Next we consider the one-dimensional gas, where the interaction is restricted to nearest 
and next-to-nearest neighbors. In this case the partition function takes the form 

N 



z^{L) = Jd^,...J deiv5(^-Ee 

n n «=1 



X 





N 

n/te)/(o + e,+i). 

As in the previous section we first compute the Laplace transform of the partition function, 

-tL 



QNit) '■= J dL Zn{L) e 


00 00 ^ 

= / d6 . . . / de^ n e-*^^/(e,)/(e, + . 

n n 7 = 1 
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(38) 



To evaluate this integral we introduce the transfer operator 

so that 5'Ar(t) can be rewritten as 
gM{t)=tiK'' 

oo oo 

^|d6.../ de;vi^(ei,6)^(6,e3)..- 



...ir(e^„i,e^)if(e^,ei). (39) 

The operator K{^, is real symmetric and therefore admits the eigenbasis expansion 

K{i,0 = Y.\,ct>M)m') (40) 

3 

with real eigenvalues \j and eigenfunctions 

00 



The eigenfunctions are normalized according to 



(41) 



and we choose the ordering of the eigenvalues by decreasing magnitude (Aq > Ai > . . . 
Consequently, in the large N limit the Laplace transform of the partition function, eq. (^ 
reduces to 



gN{t)= Ao(t) 



N 



(42) 



where we have explicitly indicated the t-dependence of the eigenvalue. Again we perform 
the Laplace inversion in saddle point approximation, which results in t being fixed to t = c. 



N 



Lc 



Zn{L) ~ [Ao(c) 
where c is determined from the saddle point equation 

AH- 1^.0. 

Ao(c) 



(43) 



(44) 



Next we calculate the joint probability distribution of n consecutive spacings, eq. (|TH). It 
turns out to be convenient to write the eigenfunctions in the form 



0(t;O = /m)e-*«/V(t;O 



(45) 
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(again the dependence on t is explicitly indicated). We express the Laplace transform 
of the integral on the right hand side of eq. ( |T8|) in terms of the transfer operator K, 
use the eigenfunction expansion eq. (^0[), perform the Laplace inversion in saddle point 
approximation (the saddle point equation is identical with eq. (0)), and arrive at 



p(6 



1 



n-l 



/(a)nm)/fe+w 



(46) 



With the help of eq. ( |46| ) we can now calculate all correlation functions. The simplest is the 
nearest-neighbor spacing distribution 



Pis) = A 



(c; sA) 



(47) 



To calculate the asymptotic behavior of the number variance we first take the Laplace 
transform of the n-th nearest-neighbor spacing distribution eq. (^O]) which yields 



'Xjjc + t/AY 
An(c) 



X 



de0o(c;O0.(c + t/A;Oe-*«/'^ 



(48) 



Hence, by eq. ([TqD , the Laplace transform of the two-level correlation function takes the 
form 



92{t) = E 



Ao(c) 



y Xo{c)~X,{c + t/A) 

oo 



(49) 



For small t this becomes 



92{t) 



1 1 A^'(c) 



t 2A X'oic) 



l + 0{t) 



(50) 



where we have used the saddle point equation (|^) and 

1 °° 

-|dee0o(c,o' = i, 



which follows from eq. ( ^Tf ) and the definition of the mean level spacing A (implying 
JdssP{s) = 1). 

We now specialize to the random matrix interaction eq. ([T^ and restrict ourselves to the 



case of P integer. The integral equation (|4T|) with the eigenf unctions in the form of eq. (^5]) 
then reads 
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J dCe-'^'e^{e + 0%it;a = A,(t)^,(t;0 • (51) 



It is clear from eq. (|5l|) that the t dependence factorizes and 



V^,(t;0 = t(^+'^/V(i;tO- (52) 



(53) 



This permits us to determine c from the saddle point equation (^41), which yields 

2/3+1 
'=^- 

Using eqs. (|23|), (0), and ( |50D we then find the asymptotic behavior of the number variance, 

1 



2(3 + 1 



In order to determine the spacing distributions we actually need to solve the integral equation 
(pT|). It is straightforward to see that it has 13 + 1 solutions each being a polynomial of degree 

^(t;0 = E%e- (54) 

The coefficients Oj are easily obtained by substituting eq. ( |5^ into eq. (|5TD. The spacing 
distributions then follow from eq. (^) and eq. (|20|) . Table || shows the explicit expression 
for the nearest and next-to-nearest-neighbor spacing distribution for /3=1, 2, and 4. 



C. General case 

Finally we extend the calculation of the previous section to an interaction between an 
arbitrary number k of neighboring particles. In this general case the Laplace transform of 
the partition function ( ]T^ ) takes the form 



N 



e X 



gM{t)= |dei.../de^ n 

i=i 

f{ii)f{i, + e,+i) ■ • • /(e, + • • • + ij+k^i) ■ (55) 

Again we seek to express gN{t) as the trace of the A^-th power of a transfer operator. To 
this end we define |jl^ 



A'(x,x') =5{x2- x[) 5{x'i - x'2) ■ ■■5{xk-i - x'^^^) x 
e-'-^l^^fJ{^)^f{x^ + X2)- ■ ■ v//(xi + ...+x,._i) X 

f{Xi + ...+ Xk-1 + 4_i) ^f{x[ + ...+X'j^_^) ■ ■ ■ X 

^fK-2 + x',^i)^m_,)e-'<-^/' , (56) 
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where x = (xi, . . . , Xk~i) and x' is defined likewise. One then easily verifies that 

= y dxi . . . y dxiv /^(xi, X2) X 

ir(x2,X3) . . .ir(xAr,Xi) . 

It is clear from the definition ( p6D that the transfer operator obeys the symmetry relation 

ir(x,x')=K(x'^x^), (57) 

where we use the notation x-^ = {xk-i, . . . ,Xi). The presence of this symmetry permits the 
eigenbasis expansion 

Ki^,^') = J2\,<P^i^)<P,i^''^), (58) 

j 

with the eigenvalues Xj and eigenfunctions 0j(x) obeying 

J dx' K(x, x') 0j(x') = Aj0j(x) (59) 

and the normalization 

J dx0j(x) (pfix^) = 5jj' . (60) 

The above considerations show that in the limit of large the Laplace transform of the 
partition function, eq. (|55D, can again be expressed as 

r n 

9N{t)= Xo{t) , (61) 



where Aq is the largest eigenvalue of the transfer operator ( p6D which parametrically depends 
on t. Hence the form of the partition function and the saddle point equation will also be 
the same as in the last section, see eqs. ( ^3]) and (|4^). 

In the following it will be useful to introduce the functions ipj{t] x) by writing 0j(t; x) in 
the form (again we indicate the t-dependence explicitly) 



where 



0,(t;x) =V^,(t;x)A/i?(t;x) , (62) 



k-l 

R{t; x) = exp ( — t Xj 

n n f{xi+...+x,+,). (63) 

j=0 i=l 



From eq. (p9|) we obtain the equation 
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/ d^k e-'^^fi^k) fi^k + e^-i) . . . fi^k + . . . + 6) X 



^j{t;^2,...,^k)=Xj{t)^,{t;^,,...,^k.,) , (64) 
which in conjunction with the normahzation condition 

oo 

J dx R{t; x) ^pj (t; x) ^pj, (t; x^) = 5jf (65) 



determines ipji^t^x.) and Xjit). 

For the calculation of the joint probability distribution of n consecutive spacings, eq. (f8|), 
we follow the same steps as in the last section, except that we have to distinguish two cases, 
namely n < k — 1 and n > k — 1. In the former case. 



oo oo 





^o(c; 6, . . . , ik-i) 0o(c; • • • , ^i) , (66) 



where c is determined from the saddle point equation (^), while in the latter case, 

n 

P(6, ■ ■ • , ^n) = Ag Vo(c; in,..., in-k+2) 6 

k-m-j . 

n n /(^- • • - ^^+.) ^o(c; 6, • • • , Cfc-i) • (67) 
i=o i=l ^ 

Eqs. ( |66|) and (|67D in conjunction with eq. (pOf ) allow the calculation of the spacing distri- 
butions. For example, the nearest-neighbor spacing distribution becomes 

oo oo 

P{s) = Jd^i.-.J da-i5(s-0/A) 



0o(c; 6, . . . , ^fc-i) 0o(c; • • • , 6) , (68) 

where j may take any value from 1 to k-1. In order to determine the asymptotic behavior of 
the number variance, we follow the same procedure as in the last section, i.e. we calculate 
the Laplace transform g{n,t) of the n-th nearest-neighbor spacing distribution, sum this 
over n to obtain the Laplace transform of the two-point correlation function, see eq. (|T9|), 
and find its small-t asymptotic, which determines the asymptotic behavior of the number 
variance, see eq. (|2^. For n < k — 2 we get 

n + l 



g(n, t) = j dil... j din+i e '=1 j d^„+2 . . . 



d^k-i 0o(c; ^1, . . . , ik-i) 0o(c; ik-i, • • • , 6) 
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whose behavior for small t is simply g{n,t) = 1 + 0{t) due to the normalization eq. (|60D. 
For n > — 2 we find 



/-^7(c + t/A)\n-fc+2r f f 



fc-i 



Summing over n and expanding asymptotically for small t we obtain 

where we have used the saddle point equation (^i]) and / ds sP{s) = 1 with given by 
eq. ( |68D . Since the k — 2 remaining terms each contribute 1 + O (t) , we finally get 

Note that this expression depends on k only by the largest eigenvalue of the corresponding 
transfer operator and it is identical with eq. (0). 

At this point we again specialize to the random matrix interaction eq. (|1^) with integer 
(3. The integral equation ( |64D which determines Xj{t) and iljj{t;x.) then reads 

oo 

/ d^k (6 + ^k^i)^ ... (6 + ... + ^i^e-'^^ X 



tPj{t- ^2, . . . , a) = A,(t) ijj{t- ^1, . . . , , (70) 
The t-dependence of the eigenvalues and eigenf unctions factorizes and 

A,(t)=r^^-%(1), 

V^,(t;x)=t('=-i)('='3+2)/4^(l;tx). (71) 
This permits us to determine c from the saddle point equation (^41), which yields 

c^^^. (72) 

Using eqs. (p3|), (p^), and (|69D we then find the asymptotic behavior of the number variance, 

S2(L)~x^, X=j^^ (L»l). (73) 



In critical statistics this quantity is related with a certain (multi)fractal dimension W7 



where d is the dimensionality of the system and rj = d — D2 and D2 is one of the multifrac- 
tional exponent defined through the behavior of the mean inverse participation ratio with 
the size of the system, L, 
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d''x\^n{x)\'^ ^ L-^^ (75) 

In our model geometrical interpretation of non-zero values of x remains unclear. 

In order to find the spacing distributions explicitly we need to determine ipQ and Aq from 
eq. ( [70D and then use eq. (p6D or ( |67D together with eq. (^). The solutions of eq. (|70D have 
the form 

P 3/3 (fc~l)fc/2 

^(t;6,---,e.-i)= E E--- E 

n=Oi2=o «fc-i=o 

Ojii2...ifc_i ^i^^2^ • • • ik-i ■ C^s) 

The coefficients aj^j2.,,jj._j must be determined numerically by substituting eq. ([76|) into 
eq. (|70D. In tables |I|, and |T1| the explicit form of -P(s) and P(l, s) are shown for fc=l, 2, 
and 3 respectively, where in each case the distributions for (3=1, 2, and 4 are given. 

These calculations become more and more tedious as or /? increase. However it is 
straightforward to find the large s asymptotic of the spacing distributions, 

P(n, s) ~ e-^^ , A = k(3 + l (s > 1) . (77) 

Note that in our model the asymptotics of the spacing distributions and the number variance 
are related as follows, 

A = -. (78) 

X 



For critical distribution of of the Anderson model at MIT point it is often assumed |18] that 



A - (79) 



To derive this relation it was assumed [113] that the probability to find n levels in an interval 



which contains in average L levels has the Gaussian form 

P„(L)~exp(-^^^), (80) 

where T!^{L) as above is the number variance. 

Therefore the probability that there is no levels inside this interval (i.e. the nearest- 
neighbor distribution) is 

Po(^)~exp(-^^^). (81) 

Assuming that, when L — > oo, T?{L) xL one gets Eq. (|79D. 

But it is clear that the assumption that Pn{L) has the Gaussian form (|80|) even for small 
n is an oversimplification and cannot be strictly valid in general. In our short-range plasma 
model Pn{L) is like the one for the Poissonian process which gives Eq. (|78|) . In the Gaudin 
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model with (3 = 2 there is no simple relation between A and x (see Eqs (|TOD and (0) (though 
for small x o^ie obtains Eq. (|79D). 

This difference between the Gaudin model when % ^ and the relation ( |78D derived 
in our short range plasma model is probably related with the different truncation of the 
interaction between two levels. In the former model the distance between two levels is 
important while in the later one only the number of levels between the two given levels is 
taken into account. 

Numerical calculations of the Anderson model at the MIT point and certain analytical 
arguments |]2B| based on the kernel ([T^) are in the favor of the relation ( |7DD at least in the 
limit of small x- The short range plasma model demonstrates that there exist (mathematical) 
models of intermediate statistics which do not obey this relation. 



III. SPECTRAL STATISTICS OF PSEUDOINTEGRABLE BILLIARDS 

Pseudointegrable systems, as introduced by Richens and Berry |T^, are dynamical sys- 
tems which are neither integrable nor chaotic. The difference between integrable and pseu- 
dointegrable systems is that for the former the phase space is foliated into 2-dimensional 
surfaces with genus g = 1 (i.e. tori) while for the latter it is foliated into surfaces of higher 
genus. The simplest example of pseudointegrable systems is plane polygonal billiards whose 
angles are all rational multiples of vr. For these models the genus of the corresponding surface 
is given by 

«-l + fE^. (82) 

k 

where mkir/rik are the interior vertex angles of the polygon, N is the least common multiple 
of the integers n^, and the sum is taken over all vertices of the polygon. The polygons with 
(7 = 1 (e.g. rectangle, equilateral triangle) are integrable, whereas those with g > 2 are 
pseudointegrable. In Refs. the spectral statistics of a number of pseudointegrable 

polygonal billiards have been investigated numerically (and in the last reference even exper- 
imentally using a microcavity) and it was found that they display level repulsion but deviate 
from random matrix theory. 



In a number of different models had been considered numerically which clearly 
demonstrate the existence of intermediate statistics. In this Section we present a more 
detailed analysis of the data. 

We consider the sequence of pseudointegrable billiards (with Dirichlet boundary condi- 
tions) in the shape of the right triangle with one angle equal vr/n, where = 5, 7, 8, 9, . . . , 30 
(the triangles with n = 3, 4, 6 are integrable) and using the boundary integral method we 
have calculated the first 20000-30000 levels for all these triangles (for some triangles even 
up to 80000 levels). 

Figure |^ shows the difference between the cumulative nearest-neighbor spacing distribu- 
tion 



N{s) = / dxP{x) (83) 

JO 

(calculated from the levels 5000-20000) for all triangles and the Semi-Poisson prediction 
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Ni{s) = l-(2s + l)exp(-2s) 



(84) 



which is obtained from our model with (3 = k = 1 (see Eq. (^)). 

Roughly four close groups of curves are observed. For the first group, consisting of the 
curves corresponding to triangles with n=5, 8, 10, 12, this difference is less than one percent 
and consequently the spectral statistics of these billiards is quite well described by the Semi- 
Poisson model. The remaining three groups correspond to the triangles with ?7,=7, 14, 18, 
with n=9, 16, 20, 24, 30, and the rest. 

The first conclusion from this and many others figures is that spectral statistics of pseu- 
dointegrable systems is not universal and depends on the billiard angles. The grouping of 
the triangles which emerges from Fig. ^ does not agree with the classification of the billiards 
according to their geometrical genus g. In first approximation the spectral statistics of the 
above triangular billiards are reasonably well described by a quantity q which we proposed 
to call 'arithmetical genus' 

\0(n)/2, neven " ^^^1 

Hence q is equal to the 'geometrical genus' g for n odd but for even n it is given by half the 
Euler function </'(n) (equal to the number of integers not exceeding and relatively prime to 
n). The first group now correspond to the triangles with q = 2, the second to g = 3, the 
third to g = 4, and the last to g > 4. It seems that for g > 4 the spacing distribution does 
not change noticeably any more, but the resulting distribution may differs from the Wigner- 
Dyson spacing distribution. We stress that this classification is only an approximative one 
and serves mostly for the crude arrangement of the spectral statistics of different triangles. 

In Figure ^ we present the evolution of the cumulative nearest-neighbor distribution for 
the triangular billiard with n = 12 with increasing energy. It is clearly seen that for higher 
energy the spacing distribution moves closer to the Semi-Poisson result (though a limiting 
distribution may deviates from the Semi-Poisson one). 

In Figure ^ the difference between the cumulative next-to-nearest distribution for the 
same triangle and the Semi-Poisson prediction (see Table D is plotted. Once more we 
observed that this distribution is close to the Semi-Poisson result and that, with increasing 
energy, the agreement is better. 

The main conclusion of this Section is that for certain pseudointegrable systems the 
short-range spectral statistics is very close to the Semi-Poisson statistics (though theoretical 
explication of this fact is not yet clear). 

IV. SUMMARY AND DISCUSSION 

In order to model intermediate spectral statistics, we considered a one- dimensional gas 
of energy levels interacting via a logarithmic pair potential, whose action is restricted to a 
small number of nearest neighbors. As shown in section 0, its correlation functions can be 
calculated analytically, so the model deserves interest as a simple reference point for compar- 
isons with numerics. In section |IT1| we performed a comparison with the spectral statistics 
of a number of pseudointegrable billiards and demonstrated that the plasma model with 
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nearest-neighbor interaction (the Semi-Poisson model) gives an excellent phenomenologi- 
cal description of the short-range spectral observables of certain pseudointegrable systems. 
Unfortunately a full theoretical understanding of the spectral statistics of pseudointegrable 
systems is still lacking. 

Of course there are many ways to interpolate between Wigner-Dyson and Poisson statis- 



/9 



tics. For the spacing distribution a well-known interpolation is the Brody distribution p2 
which shows fractional level repulsion. Since for intermediate spectral statistics P{s) ~ s 
for small spacings, the Brody distribution cannot give an adequate description. A very 
natural way to construct interpolating ensembles is to take a weighted average between the 
ensemble of diagonal random matrices and one of the standard random matrix ensembles, 
e.g. for /5 = 1 

H = -f^Poisson + A i^GOE (86) 



(see e.g. ref. ||2^). However, as we checked numerically, for no value of A the resulting 
distribution is close to the distribution (|^), which we have found to be an excellent fit to the 
spacing distribution of the pseudointegrable triangles with q = 2 and also to that of several 



other systems []TU|. Hence this interpolation do not seem to be suitable for the description 
of intermediate statistics. 

All data suggest that a necessary requirement of theoretical description of intermediate 
statistics is the screening of two-body potential. Its precise form is not yet known (but see 



24|- |26]) and may depend on the problem considered. 
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FIGURES 




FIG. 1. N particles on a circle of circumference L. The positive spacings between consecutive 
particles measured along the circle are denoted by with j = 1 . . . N. 
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FIG. 2. Difference between the cumulative spacing distribution for the pscudointegrable ratio- 
nal right triangles with n < 30 (calculated from the levels 5000-20000) and the Semi-Poisson curve. 
See explanation in the text. 
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FIG. 3. The same as in Fig. |2| but for n = 12. The Ihin line corresponds to 1-34000 levels and 
the thick one includes 34001-68000 levels. 
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FIG. 4. The difference between the cumulative next-to-nearest spacing distribution for n = 12 
and the Semi-Poisson prediction. Two curves are the same as in Fig. ^ The dotted line corresponds 
to GOE 
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TABLES 



TABLE L The spacing distributions P{s) and P(l,s) for the one-dimensional gas model with 
nearest-neighbor interaction. 

{k = 1) (3 = 1 p = 2 (3 = A 

P{s) 4se-^' 27/2 s^e-^'^ 3125/24 s^e'^^ 

P{l,s) 8/3 s^e-'^' 243/40 1953125/72576 s^e"^* 



TABLE II. Same as Tab. |, but with nearest- and next-to-nearest-neighbor interaction. 



{k = 2) 



f3=l 



(3 = 2 



P(l, S) = s3/3+le-{2/3+l). ^ . ^.^j 



Co = 


2.4773 


Co 


= 6.2603 


Cl = 


6.0681 


Cl 


= 24.9958 


C2 = 


3.7159 


C2 


= 41.1181 






C3 


= 32.2768 






C4 -■ 


= 10.4386 


do = 


2.5054 


do 


= 5.9771 


di = 


: 3.068 


dl : 


= 11.9326 


d2 = 


0.7516 


d2 


= 9.5151 






ds 


= 3.3018 






di 


= 0.4746 





Co = 35.6018, Cl = 252.5852 
C2 = 826.3343 , cg = 1630.9839 
C4 = 2127.3725, cg = 1880.5911 
C6 = 1103.6323, C7 = 395.1819 

C8 = 66.6078 
do = 29.4495 , di = 104.4681 
d2 = 169.5742, dj, = 164.1739 
di = 103.3668 , 4 = 43.2920 
= 11.9204, dj = 2.0002 
ds = 0.1586 



23 



TABLE in. Same as Tab. |, but with three interacting neighbors. 

{k = 3) 13 = 1 13 = 2 /? = 4 

Co = 2.1342 Co = 4.9711 co = 24.0711 , ci = 282.9984 

P[s) = s/3e-(3/^+i)'* Y.j CjS^ ci = 7.8196 a = 31.6156 C2 = 1615.27 , C3 = 5952.78 

C2 = 11.6945 C2 = 93.9581 C4 = 15890.11 , C5 = 32666.46 

C3 = 9.0743 C3 = 171.4838 ce = 53678.64 , C7 = 72222.20 

C4 = 3.6855 C4 = 212.8904 cs = 80836.91 , cg = 76044.48 

C5 = 0.6259 C5 = 188.1012 cio = 60480.73 , cn = 40759.87 

C6 = 120.0046 C12 = 23247.21 , C13 = 11165.36 

C7 = 54.6739 C14 = 4473.73 , C15 = 1472.98 

cs = 17.0418 C16 = 389.3105 , C17 = 79.6241 

eg = 3.2829 ci8 = 11.8577, cig = 1.1466 

cio = 0.2968 C20 = 0.0541 

do = 1.1340 do = 1.6110 do = 2.9196 , di = 29.5956 

P(l, s) = s3/3+ie-{3/3+i)s ^.gj ^ 3 g334 ^ 8 3782 ^2 = 146.1401 , ^3 = 467.9231 

d2 = 5.0189 d2 = 23.0847 ^4 = 1090.9150 , = 1971.2442 

d3 = 3.7112 ds = 37.4284 de = 2868.9524 , dr = 3449.8348 

d4 = 1.5037 d4 = 42.0839 4 = 3488.0347 , dg = 3001.6809 

4 = 0.3100 (is = 34.4593 dw = 2216.98 , du = 1412.97 

= 0.0259 de = 20.9804 di2 = 779.3162 , d^ = 372.3224 

dj = 9.5327 di4 = 153.8623 , dis = 54.8050 

ds = 3.2010 di6 = 16.7526 , du = 4.3271 

dg = 0.7747 di8 = 0.94224 , dig = 0.1687 

dio = 0.1284 d2o = 0.0241 , d2i = 2.76 * lO'^ 

dii = 0.0131 d22 = 2.5 * 10"^ , d23 = 1.9 * lO"'' 

di2 = 6.2 * 10~^ d24 = 6.2 * lO"'^ 
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